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Abstract 

We study control problems for linear systems in the behavioral 
framework. Our focus is a class of regular controllers that are equiva- 
lent to the canonical controller. The canonical controller is a particular 
controller that is guaranteed to be a solution whenever a solution exists. 
However, it has been shown that in most cases, the canonical controller 
is not regular. The main result of the paper is a parametrization of 
all regular controllers that are equivalent to the canonical controller. 
The parametrization is then used to solve two control problems. The 
first problem is related to designing a regular controller that uses as 
few control channels as possible. The second problem is to design a 
regular controller that satisfies a predefined input-output partitioning 
constraint. In both problems, based on the parametrization, we present 
algorithms that does the controller design. 

1 Introduction 

In this paper, we discuss control problems for linear differential systems in 
the behavioral approach. The behavior of a system is the set of trajectories 
that are compatible with the laws that describe the system. In the continu- 
ous time case, the behavior is the set of solutions of the differential equations 
that describe the system. In the discrete time case, it is the set of solutions 
of some difference equations. 
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Figure 1: Control in the behavioral approach. 



Standard control problems in the behavioral approach to systems theory 
can be formulated as follows [HE1E]. Given is a plant to be controlled that 
has two kinds of variables: to-be- controlled variables and control variables. 
A controller is a device that is attached to (or an algorithm that acts on) 
the control variables and restricts their behavior. This restriction is im- 
posed on the plant via the control variables, such that it (indirectly) affects 
the behavior of the to-be-controlled variables (see Figured]). The resulting 
behavior is called the controlled system. 

As part of the control problem, one is given a specification, which is 
expressed in terms of the to-be-controlled variables. The objective of the 
control problem is to make the controlled system satisfy the specification. 
If there exists a controller such that this objective is satisfied, we say that 
the specification is implementable. 

In [U [5], a particular controller design, called the canonical controller 
was introduced. This design has the nice property that it implements the 
desired specification if and only if the specification is implementable. How- 
ever, analysis on the regularity of the canonical controller reveals that it 
is maximally irregular [6]. Regularity is a desirable property for the in- 
terconnection [TJ [2], which we will explain in Section [2j In this paper, we 
show that there exist regular controllers that are equivalent to the canonical 
controller, and we provide a parametrization of all such controllers. This 
parametrization is then used to solve two control problems: 

1. The problem of control with minimal interaction [7]. This prob- 
lem is about designing a regular controller that interacts with the plant 
with as few control variables as possible. The motivation behind this 
problem is as follows. Consider a situation where the plant and the 
controller are separated by a large physical distance. We need a com- 
munication link between the plant and the controller to establish the 
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interconnection. It is therefore favorable to have as few control vari- 
ables as possible, so that the amount of communication links/channels 
can be minimized. 

2. The problem of control with I/O partitioning constraint. This 
problem is about designing a regular controller that respects a con- 
straint on the a priori partitioning of the control variables into input 
and output variables. 

The results in this paper are presented in the form of continuous time 
systems. However, they also hold for discrete time systems, as we replace 
the differential operator with the discrete time lag operator. 

2 Background material 

For linear differential systems, the plant is typically described as a set of 
linear differential equations that relate the variables. Throughout this paper, 
we denote the control variables as c and the to-be-controlled variables as w. 
The dimensions of c and w are denoted as c and w respectively. A behavioral 
model of the plant system that captures the relevant relation between w 
and c is called the full plant behavior, and is denoted by Pfuii- The full plant 
behavior can be compactly represented as follows. 



where R and M are polynomial matrices with appropriate dimensions. We 
denote the class of polynomial matrices with indeterminate £, g rows, and 
q columns over the real field as R ffX,? [£]. 

The representation of the behavior in the form of (pQ) is called a kernel 
representation, the reason being that the behavior is simply the kernel of a 
linear differential operator. Kernel representations of a given behavior are 
not unique. For example, if U(d/dt) is a linear differential operator, whose 
kernel consists of only the zero trajectory, then the kernel of U{d/dt) o 
R(d/dt) is the same as that of R(d/dt). Square polynomial matrices U(£) 
such that 



are called unimodular matrices. It can be proven that the inverse of U(£) is 
also a polynomial matrix. A stronger result that relates unimodular matri- 
ces and equivalent kernel representations is that two kernel representations 






(2) 
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R\(d/d£) and R2(d/dt) with the same number of rows are equivalent if and 
only if there is a unimodular matrix U(£) such that -Ri(C) = ^(0^2(0- 

Although the kernel representation of a behavior 03 is not unique, there 
is a unique integer p(03), which is the minimum number of rows a kernel 
representation of 03 can have. This number is also the row rank of any 
kernel representation of the behavior. A kernel representation with the 
minimum number of rows (i.e. equal to its row rank) is called a minimal 
kernel representation. The number p(03) is called the number of outputs of 
23. 

Suppose that a behavior 03 is given by 

»:={.|*(|)» = o}, (3) 

where R is full row rank and has p(03) rows. We can partition the variables 
in w into wi and W2 such that ([3j) becomes 

<8 := Uwuim) \R 1 (J^w 1 + R 2 (J^ W2 = o\, (4) 

where R\ is a square full row rank polynomial matrix. Such a partition is 
called an input-output partition where wi is the output and W2 is the input 
to the system. Notice that the number of outputs of 03 is p(03). 

In this paper, we restrict our attention to infinitely differentiable func- 
tions. Thus, the full plant behavior consists of all signal pairs (w, c) that 
are strong solutions to the kernel representation ([1]) |8j . 

V iull :=^w,c)G^(R,W +c ) \r(J^ W + m(J^Jc=o\. (5) 

If we eliminate the control variables from the full behavior, we obtain the 
so called manifest behavior, which is denoted by V . Thus, 

V:={we £°°(IR,IR W ) | 3 c G £°°(M,R C ) such that (w,c) G V m }. (6) 

If we rewrite the kernel representation ([!]) as 



Mi) M i(i) 





w 




c 



0, (7) 



where M\ and R 2 are full row rank matrices, then the manifest behavior V 

'A' 



is the kernel of R 2 (4) (cf. [8] Chapter 6). 
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A controller C is a behavior containing all signals c allowed by the con- 
troller: 

C:= jcee°°(M,R c ) | C (Jl) c = \- ( 8 ) 
The controlled behavior is then defined as 

K:={w G £°°(M,R W ) | 3 c G £°°(IR,M C ) such that (u>,c) G P M1 and c G C}. 

(9) 

The controlled behavior K, is obtained by eliminating the control variables 
from the following kernel representation. 

*(D"' +M (!H ' 

c(^-)c = 0. (10) 



The specification S is given by the following kernel representation 

s (I).- a (ii) 

The objective of the control problem is to find a controller C such that 
KL = S. If such controller exists, then S is said to be implementable and the 
controller C is said to implement S. 

Clearly, the implementability of a specification S is a property that de- 
pends on the specification itself as well as the plant. The following result is 
proven in [9j [10] . 

Theorem 1 (Willems' lemma) Given Vf u u as a kernel representation of 
|T]). A specification S is implementable if and only if 

NQSQV, (12) 

where M G £ w is the hidden behavior defined by 

N := {w G e°°(R,lR w ) | (w,0) G V full }. 

Quite often, in addition to requiring that the controller implements the 
desired specification, we also require that the controller possesses a certain 
property with respect to the plant. A property that has been quite exten- 
sively studied is the so called regularity [HI El HI Q3] . A controller 

cer(l,l c ) I c(^-)c = o\, (13) 
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where C is full row rank, to be regular if 



rank 



R M 
C 



rank [ R M ]+ rank C. 



(14) 



It can be proven that nonregular interconnections affect the autonomous 
part of the systems [I], which, in many cases would be undesirable or unre- 
alistic. 

Remark 2 Although the characterization of regular controllers suggests that 
regularity is a representation dependent property, it is actually not. Notice 
that Ji^D is equivalent to saying that the number of outputs of the controlled 
system is the the sum of those of the plant and the controller. The number of 
outputs of a system, as discussed earlier in this section, is a representation 
independent quantity. The interested readers are referred to JU \12\ \13\j for 
more discussion on the behavioral interpretation of regularity. 

If the specification S is such that there exists a regular controller C that 
implements it, then S is said to be regularly implementable. Necessary and 
sufficient conditions for regular implementability were derived in [2]: 

Theorem 3 Given the full plant behavior Vfull- A specification S is regularly 
implementable if and only if 

1) it is implementable, i.e., M C S QV and 

2) S + V ctr = V. 

The symbol V ctr denotes the controllable part of the manifest behavior V. 

3 The canonical controller and its regular equiva- 



In this section, we review the idea of canonical controller and its properties. 
Given a full plant behavior V{ n \\ and a specification S, the behavior of the 
canonical controller C can is defined as 



C can := {c G R c ) | 3 w € £°°(M, IT) such that (to, c) G P M and w G S}. 



A kernel representation of the canonical controller can be obtained by elim- 
inating w from the following kernel representation 



lences 



(15) 




(16) 



The canonical controller has the following property. 
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Theorem 4 (cf. f5jj) The canonical controller C can implements the specifi- 
cation S if and only if S is implementable. 

We define the control manifest behavior of the plant, V c as 



V c := {c e £°°(R,M C ) \ 3w<E (£°°(M, M w ) such that (w,c) G Vmi}- (17) 



A kernel representation of V c can be obtained by eliminating w from the 
kernel representation of Pfuii- The canonical controller has the property of 
being least restrictive in the following sense. 

Proposition 5 (cf. Assume that the specification S is implementable. 
For any controller C that implements S, we have that 



Thus, any trajectory of the control variables of the plant allowed by C is also 
allowed by C can . 

Another important property of the canonical controller that is relevant to 
our discussion in this paper, is that it is maximally irregular, in the following 
sense. 

Theorem 6 (cf. f^) Assume that the specification S is implementable. 
The canonical controller C can is regular if and only if every controller that 
implements S is regular. 

Although the canonical controller is maximally irregular, there are regu- 
lar controllers that are equivalent to it. By equivalent controllers, we mean 
the controllers that allow the same set of c trajectories of the plant as the 
canonical controller does. The class of such controllers is defined as follows. 

Definition 7 The class of regular controllers that are equivalent to the 
canonical controller is denoted as Ccan, and is defined as 



The following theorem provides a necessary and sufficient conditions for 
the nonemptyness of the class £ C an- 

Theorem 8 The class is nonempty if and only if the specification S is 
regularly implementable. 



(C n v c ) c (C can n v c ) ■ 



(18) 



CI -={C\C is regular and (C n V c ) = (C can n V c )}. 



(19) 
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Proof. The (only if) part of the theorem is obvious. We shall prove the 
(if) part. Suppose that S is regularly implementable. There exists a regular 
controller that implements S. We denote this controller as C. By definition, 
we have that 

(a) For all w€5, there exists a c G C such that (w, c) G Pfuii- 

(b) For all c G C, (w, c) G Pfuii implies u> G 5. 

Define another controller 

C':=C + C can . (20) 

We shall prove that C G £ C an, that is 
(a') C' is regular. 

(b') Cnv c = c can n P c . 

The statement (a') follows from the fact that C C C and the regularity 
of C. To prove (b'), first we show that C implements S. From here, (b') 
follows from the fact that C can C C and the property of C can being the least 
restrictive controller (see Proposition [5]) . 

Showing that C implements S means showing that 
(a") For all w G 5, there exists a d G C such that (w, c') G Pf u u. 
(b") For all c' G C, (w, c') G V in \\ implies w G S. 

Statement (a") follows from (a) and the fact that C C C. To show that 
(b") holds, notice that any c' G C can be written as c + c can with c G C and 
Ccan G C can . Also notice that for all c can G C can , there exists a ui can G 5 such 
that (-u; can ,c can ) G Pftji- Thus, 

(w, c') G Pmi =^ (w - w can + w 

can; <- ~r i-can ) G P fu ii 

linearity . . 

=> ((lO - U> C an),c) G Pfull 

(w - u> can ) G 5 

linearity _ ^, 
=> U> G O. 



The proof of Theorem [8] also implies the following important property 
of <f reg 

ui ^can- 

Theorem 9 Given a control problem with a regularly implementable speci- 
fication S. If C is a regular controller that implements S, then there exists 
a regular controller C G Ccan that implements S and C C C' . 
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One of the main results presented in this paper is the parametrization of 
all controllers in £can- Before we can obtain the parametrization, we need 
the following lemma. 

Lemma 10 Let a plant V be given as the kernel of a full row rank R (^) 
and a regular controller C be given as the kernel of a full row rank C (^) . 
Denote the full interconnection 

JC:=VnC. 

Let (L/c denote the set of all controllers (not necessarily regular ones) that 

(i) have at most as many outputs as C and 

(ii) also implement K, when interconnected with V . 

A controller C £ if and only if its kernel representation can be written 
as VR + C for some matrix V. Moreover, every controller in C £ £)c has 
the following properties. 

(a) C is regular. 

(b) C has exactly as many outputs as C. 

Proof, (if) Suppose that a controller C is the kernel of {VR + C), then 
V fl C is given by the kernel of 



(21) 



This shows that VnC' = VnC = IC. Moreover, since C is a regular con- 
troller, it follows that (VR + C) is a minimal kernel representation of C . 
Thus, properties (a) and (b) are verified. 

(only if) Suppose that a controller C satisfies (i) and (ii) above. This 
controller can be written as the kernel of a matrix (not necessarily minimal) 
C'(^) with as many rows as C(^). We know that there is a unimodular 
matrix U such that 



R 




' L " 




' R ' 


VR + C 




V L 




C 



' R ' 




' U u 


C/12 " 




' R ' 




' R ' 


C 




_ U 2 i 


U22 _ 




C 




a 



u 

We shall prove that we can assume U to be of the form 

U 



L 
V L 



First, we find a unimodular matrix W such that 

RW =[ D ] , 



(22) 



(23) 



(24) 
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where D is a square nonsingular matrix. We then use the following notation 



R 
C 

R 

a 



It follows that (1221) can be rewritten as 



W =: 


D 





. C i 


c 2 


W =: 


D 





. c[ 


C'2 . 



(25) 
(26) 



U 



D 

C\ C*2 



and since W is unimodular, 
Consequently, we have the following equations 



D 









c 2 





^21 ^22 










w 


-1 






d 


G 2 








D 







D 







Ci 


c 2 




. ci 


^2 . 



U U D + C/ 12 Ci = D, 
U l2 C 2 = 0, 

c/ 2 i£> + u 22 d = c[, 

U 22 C 2 = Co. 



(27) 



(28) 



(29a) 
(29b) 
(29c) 
(29d) 



Since the controller C is regular, C 2 must be full row r ank. Now, (j29bjl 
implies that U\ 2 is a left annihilator of C 2 . Consequently 



U 



12 



Substituting this to ()29a|) yields 



U 



ii 



0. 



(30) 



(31) 



Since U is unimodular, this implies that U 22 is unimodular. Thus, we can 
conclude that 

/ 

U 2 i U 22 



u 



(32) 



with U 22 unimodular. Furthermore, C" := U 22 C is also a kernel represen- 
tation of C so we can assume U 22 to be the identity matrix without any loss 
of generality. ■ 

Now we can parametrize the elements of £cai as follows. 
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Theorem 11 Let the control manifest behavior of the plant V c be the kernel 
of P c (jA and a controller C € £ C an be the kernel of C ( Jj) . Assume that 
both P c and C are full row rank. A controller C is also an element of £ C an 
if and only if it is the kernel of V ( Jj) P c (^) + C (^) for some polynomial 
matrix V(£). 



Proof. The full plant behavior can be represented by 



Hi) mU) 





w 




c 



(33) 



where R is full row rank. It follows that a controller C represented as the 
kernel of C ( Jj) is regular if and only if 



rank 



Pc 

a 



rank P c + rank C . 



(34) 



This is equivalent to saying that the interconnection of V c and C is regular. 
Therefore, we can apply Lemma [10] (by replacing K, with C can and V with 
V c ) and obtain the parametrization of all elements in <t T c 



r reg 
-can • 



4 Control with minimal interaction 
4.1 Problem formulation 

Consider the following definition of irrelevant variables. 

Definition 12 Let a behavior 53 be given by the kernel representation 

fii (£) roi+i K£) ro2=o - <35) 

Lf Ri is the zero matrix, then the variables in wi are said to be irrelevant 
to 53. 

Notice that whether or not some variables are irrelevant to a behavior 
is not a matter of representation. Rather, it is a property of the behavior. 
It means for every (wi,-^) G 23 we can always replace w\ by any infinitely 
differentiable trajectory w[ and have that {w' 1 ,W2) is still an element of 
23. Hence, although wi is explicitly present in the description of 23, the 
information about its trajectory is irrelevant. 
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Figure 2: Control with irrelevant control variables. 



The problem of control with minimal interaction that we are addressing 
in this paper can be formulated as follows. 

Control with minimal interaction. Given are the full plant behavior 
"Pfuii © and specification S. We assume that the specification S is regu- 
larly implementable. Construct a regular controller C that implements S 
with as many irrelevant variables as possible, or equivalently as few relevant 
variables as possible. 

The controller to be designed is called the controller with minimal in- 
teraction. When some control variables are irrelevant to the controller, we 
can realize the controller without using these variables. See Figure [2] for an 
illustration. 

4.2 The solution 

We are going to use the parametrization of £can that we derived in the 
previous section to solve the problem of control with minimal interaction. 
First, consider the following lemma. 

Lemma 13 Let a behavior 03 be given by the kernel representation 



If wi is irrelevant to 03, then it is also irrelevant to any 33' 3 23. 
Proof. The kernel representation of any 23' 5 23 can be written as 




(36) 



FR 1 — )wi + FR 2 — )w 2 = 





12 



for some polynomial matrix F. Clearly FR\ = 0, thus wi is irrelevant to 

Lemma [13] and Theorem [9] tell us that it is sufficient to search for the 
controller with minimal interaction in (£can, instead of in the set of all regular 
controllers. This is an advantage, since we can parametrize all the controllers 
in £can, as shown in Theorem [TT1 To solve the problem of control with 
minimal interaction, we need to find an element of £can with as many zero 
columns as possible. Generally, since there are finitely many columns, there 
is a maximal number of zero columns that can be attained. However, there 
is no guarantee that this number is attained by a unique controller. In fact, 
generally speaking, it is not. 

The procedure to compute a regular controller that implements S and 
has as many irrelevant variables as possible can be summarized as follows. 

Step 1. Construct the canonical controller C can for the problem. Since S 
is regularly implementable, we know that the canonical controller im- 
plements <S. 

Step 2. Construct a controller C € Ccan- The proof of Theorem [8] describes 
how to construct C from a regular controller. Denote the kernel rep- 
resentation of C and the control manifest behavior, V c , by C(-jr) and 
P{4f) respectively. 

Step 3. The kernel representation of the controller with minimal interac- 
tion can be found by finding a matrix V such that C + VP has as 
many zero columns as possible. 

The algebraic problem related to the third step has a combinatorial 
aspect in it, as we generally need to search for the answer by trying all 
possible subsets of the columns. This situation gives rise to a computational 
challenge, namely to design an algorithm that can handle this combinatorial 
problem efficiently. Before we proceed to discuss the algorithm, we establish 
an upper bound for the number of irrelevant variables that can be attained 
in the controller with minimal interaction. 

Lemma 14 The controller with minimal interaction can have at most c — 
p(C) irrelevant variables. Here c denotes the number of all control variables 
(the cardinality of c) andp(C) denotes the number of output variables inC, 
which is any regular controller that implements S. 

Proof. From the definition of regularity, we know that all regular controllers 
that implement S have the same number of outputs, i.e., p(C). This is the 
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number of rows in a minimal kernel representation of the controller. It is 
easily seen that the number of columns is c. If a regular controller has more 
than c — p(C) irrelevant variables, then the nonzero entries of its kernel 
representation form a tall matrbi0, and thus cannot be minimal. ■ 

Notation 15 In the subsequent discussion, we denote the entry on the i- 
th row, j-th column of C as C^. The j-th column of C is denoted as C 9 j 
and the i-th row as C^. In a similar fashion, we also define Vij, V,j, Vi,, 
Pij , P,j , and . Moreover, we denote the greatest common divisor of the 
polynomials in P,j as ttj . 

Notice that the j-th column of C+ VP is zero if and only if C,j + VP,j = 
0. Consider the following proposition. 

Proposition 16 There exists a V such that 

Cj + VP.j = (37) 

if and only if ttj divides C,j . 

Proof, (if) Suppose that Cij = Kiirj, where polynomial. Since ttj is 

the greatest common divisor of the polynomials in P 9 j, there exists a row 
vector v such that the Bezout identity 

V ■ P.j = TTj 

is satisfied. It follows that choosing V such that Vij = —KiV will give us 

dSZJ). 

(only if) Suppose that C,j = —VP,j. It means Cij = —Vi,P,j. Since irj 
divides P,j, it also divides Cij. ■ 

The process of constructing the matrix V that corresponds to a controller 
with minimal interaction can be considered as a recursive process. 

The recursive computation for V. Suppose that we are given C and P, 
and we want to construct V such that C + VP has as many zero columns as 
possible. Suppose that first we want to nullify the i-th. column of C + VP. 
We proceed with the following steps: 

1. Check if C., is divisible by TTj. If not, the procedure stops here, 
otherwise, we denote C.j as [ Hi^u Hi^i • • • ^i^p(c)i ] • 

A tall matrix is a matrix, in which there are more nonzero rows than there are columns. 
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2. Compute a unimodular matrix U such that P := UP is such that its 
i-th column is [ tti • • • ] . 

3. Define V := Vl/ -1 . We then have that 

C + VP = C + VP. 

It follows that the i-th column of C + yP is zero if and only if the first 
column of V is [ -k u -K 2i ■ ■ ■ -«p(C)i ] • 

4. We can write V and P as 

V=[V 1 V 2 ] , P = 

where Pi is the first row of P and V\ is the first column of V, which is now 
known. 

5. Notice that 

C + VP = C + V 1 P l + V 2 P2- (38) 

Define C := C + V\P\. The i-th column of C is zero by the construction in 
the previous steps. If we want to proceed by, say, nullifying the j-th column 
of C + VP, then the problem is to design V 2 such that the j-th column of 
C + V2P2 is zero. So now we arrived at a problem similar to the one we 
started with. However, now we have (at least) one less column to nullify and 
one less column of V to design (since the first column is determined). We 
define an algorithmic function that does the computation described above. 
Algorithm of the function (C, P 2 , V\, U, f ail, skip) = nullif y(C, P, i) 

1. Check if C.j is zero. If yes, return (C, P, 0, 1, false, true). Otherwise, 
go to step 2. 

2. Compute a unimodular matrix U such that P := UP is such that its 
i-th column is [ 7Tj • • • ] . 

3. Check if 7Tj divides C.j. If not, fail = true, return (0, 0, 0, U, f ail). 
Otherwise fail = false and go to step 4. 

4. Compute V\ = -^-C 9i . 

, where P\ is the first row of P. 
6. Compute C = C + V 1 P 1 . 




P 

5. Partition P into ~ 
P2 
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7. Return (C, P 2 , Vi, U, fail, false). 



The following example illustrates the algorithm. 



Example 17 Lei 

^(0 = 



£ £ 1 
e + i e o 



Suppose that we want to design V(£) smc/i i/iai t/ie /irsi column of C + 
is zero, i.e. i = 1. M^e ster£ itrai/t step i o/ £/ie algorithm above, and since 
the first column of C is not zero, skip=false and we go to step 2, where 
we obtain 



-1 1 



u(0 = 

HO ■= u(OP(0 



10 -1 

o -c -e-i 



Notice that the greatest common divisor of P,i is 1, which divides C,i. This 
means fail = false and we go to step 4, 5 and 6. 

m) = C-f, Pi(0 = [ io-i], 

p 2 (o = [o -e-i], <?(*) = [ o e-i e 2 -e-i] 

T7ms, w;e can verify that the first column of C has been nullified. 

The decision on the order of the columns that we nullify involves a combi- 
natorial search. Now, we are going to develop a depth-first search algorithm!! 
(see for example, |14j ) that can compute the controller with minimal inter- 
action. First, we are going to define an efficient data representation for the 
depth-first search. Consider the set A := {1,2,3, ••• , c}, where c is the 
cardinality of c, which is also the number of columns in the representation 
of the controller. We define S to be the set of increasing strings formed by 
the elements of A nonrepeatingly, such that if s £ S and |s| is the length of 
s then 

Si < Sj+i, 1 < i < \s\ . (39) 



2 Depth-first search is a standard term in computer science. It is an algorithm for 
traversing or searching a tree, tree structure, or graph. Intuitively, in this algorithm, one 
starts at the root (selecting some node as the root in the graph case) and explores as far 
as possible along each branch before backtracking. 
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Figure 3: The search tree for c = 3. 



The empty string is denoted by e and its length is defined to be 0. Thus, 
the numbers in a string s G S in increasing. There are 2 C elements of S, 
each of which represents an element of the power set of A. We can assign 
lexicographic ordering (see for example, [13]) on the elements of S and sort 
them. That is, for any two distinct elements s, s' € S, 

c 

s<s ' J2(s'i - Si)N { > 0, (40) 
i=l 

where N is any integer larger than c, and under the convention that 

Si = 0,i > \s\ . (41) 

We then define a subset T C S, by 

T := {s € S | a w = c}. (42) 

That is, T is the set of strings that end with c. For example, the elements 
of S and T for c = 3, in ascending order, are 

S = {£,1,12,123,13,2,23,3} , 
T = {123,13,23,3}. 

The search tree for the problem, where c = 3 is shown in Figure [3j Notice 
that each element of S represents a node on this tree. The ordering of S 
tells us the order in which the nodes are visited in the depth-first search. 
The elements of T represent the terminal nodes, each which represents a 
path from the initial node (the top of the graph) to a terminal node. 

We define the following operations on S. The prefix operator pre : 
S — > 2 s is such that pre(s) is the set containing all the prefixes of s. The 
operator » + and »~ are such that s + is the last symbol in the string s and 
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s is the string formed by removing the last symbol from s. The operator 
|-, •) : S x N — > T is defined as follows. 

\s, k) := the smallest element of {a G T | a > s, |er| > k}. (43) 

If the set in (|43|) is empty, then \s, k) := e. The operation |-, • j : S x S — > S 
is defined as follows. 

|s,s'| := the smallest element of pre(s') Dpre(s)^. (44) 

The operation [•] : S — > S is defined as follows. 

[~s] := the smallest s' G S s.t. s ^ pre(s'). (45) 

Denote the cardinality of S by |S| = 2 C — 1, and the i-th element of S by 
S(<). 

The following algorithm takes polynomial matrices C and P as inputs 
and returns a matrix V such that C + VP has as many zero columns as 
possible. 

Algorithm of the function V = computeV(C, P) 

1. If C already has zero columns, compute a unimodular permutation W\ 
such that C := CW\ has all the zero columns on the left. Otherwise, 
Wi = I. Define P := PW X . 

2. Denote the maximum possible number of zero columns (see Lemma 
[T4"|) as M. If the number of zero columns of C is less than M, go to 
step 3, otherwise return V = 0. 

3. Initialize the variable k = 0, n max = 0, and the strings a and <T max are 
both empty. 

4. Define C[0] = C, P[0] = P, f ail(0) = false. 

5. If n max < M and k < |S|, go to step 6. Otherwise, go to step 13. 

6. Define a = S(k). 

7. If f ail(/c) = true then change k such that S(k) = \a~\ and go to step 
9, otherwise go to step 8. 

8. If {s G T | s > a, \ s\ > n max } = go to step 13, otherwise change k 
such that S(k) = |<r, |cr, n max }|. 



18 



9. Let j be such that S(j) = (S(k))~, compute (C[k],P 2 [k], Vi[k],U[k], 
fail[fc],skip[A:]) =nullify(C[j],P[j],(S(k))+). 

10. If fail [A;] = false, then go to step 11, otherwise go to step 12. 

11. If |S(/c)| > n max then modify n max = |S(A;)| and <7 max = S(A;). 

12. Go to step 5. 

13. Initialize a = cx max . If <r max = e, then return V = 0. Otherwise go to 
step 14. 

14. If \a\ = go to step 19, otherwise go to step 15. 

15. Let k be such that S(k) = a. If skip[/c] = true then go to step 18. 
Otherwise go to step 16. 

16. V = Q law dim(C)xcoldim(P 2 [k]y 

17. V =[ Vi[k] V ] , redefine V = VU[k}. 

18. a = (S(fc)) - and go to step 14. 

19. Return V. 

Example 18 Consider the matrices given in Example \ 1 1\ 

If we apply the algorithm above to this example, then the following steps are 
going to executed. 

1. The first column of C will be nullified, as shown in Example \ 1 1\ Thus 
the maximum number of column that can be nullified by the algorithm so far 
is 1. 

2. See Figure O The algorithm is now at the first branch from the left, in 
the figure. The algorithm will subsequently try to nullify the second column. 
That is, it will try to find polynomial matrix V2 such that the second column 
of 

[0 £ - 1 e - £ - 1 ] + V 2 [ -£ -£ - 1 ] 

is zero. Since this is not possible, it will try to nullify the third column, 
which corresponds to the second branch from the left in Figure [3 This is 
also not possible. 



p(0 
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3. The algorithm then tries to nullify the second column of C (third branch 
from the left in Figure^. This is not possible since the greatest common 
divisor of P.2 is £, and it does not divide C,2 . 

4- Since the remaining branch in Figure^ consists of only one element, it is 
not possible to find a combination of columns, consisting of more than one 
column, that can be nullified. The algorithm then terminates and the final 
result is the nullification of the first column of C. 



5 Control problem with input-output partition con- 
straint 

5.1 Problem formulation 

One of the features of the behavioral approach to systems theory is that no 
a priori distinction is made between input and output variables of a system 
\ib\ [8]. This means that given a certain law that describes the system, the 
system is identified by the collection of its trajectories as is. Therefore, it 
is not necessary to have any input-output structure when describing the 
system. 

However, when two systems are interconnected, sometimes some input- 
output structure can emerge naturally as a constraint. Consider the follow- 
ing example. 

Example 19 Consider a tank filled with water as shown in Figure [7} On 
top of the tank is an inlet from which a variable flow of water can get into 
the tank. We denote the water flow from this inlet as e. On the bottom of 
the tank, there is an opening connected to a pump that can pump water out 
of /into the tank. We denote the amount of water flow pumped out of the 
tank as u. The tank is also equipped with a sensor that measures the change 
of volume of water inside the tank, the measurement of the sensor is denoted 
as d. The mathematical model of this system can be simply written as 

d(t) = e(t)-u(t). (46) 

Now consider the following control problem. Given d and u as control vari- 
ables, we want to design a controller such that the level of water is constant, 
i.e. e(t) = u(t). In other words, we aim at perfect tracking of e by u. Intu- 
itively, we know that such task cannot be accomplished. However, consider 
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Figure 4: The water tank system in Example [T9J. 



the following construction. First we write the plant behavior in a kernel 
representation. 



V = {(e, u, d) | e(t) - u(t) - d(t) = 0}. 
We then take a candidate controller C expressed by 

C = {( u ,d) \ d{t) = 0}. 
The interconnection V || C is represented by the 



1 -1 -1 
1 



e 
u 
d 



0. 



(47) 



(48) 



(49) 



Notice that the interconnection exhibits the following features. 

(i) The interconnection is a regular interconnection. In fact, it is even a 
regular feedback interconnection fj]/. 

(ii) The controller is expressed only in terms of u and d. 

(Hi) In the controlled behavior, perfect tracking e(t) = u(t) is attained. 
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Figure 5: The water tank system in Example [20l 



In the example above, the proposed controller is regular and accom- 
plishes the control task. However, this is still counter intuitive, and impos- 
sible to implement. The variable d is a measurement coming from a sensor, 
and yet we use it to enforce control on the system. Otherwise stated, we 
control the system by restricting the reading of a sensor. Now, consider the 
following modification of the example. 

Example 20 Let us swap the name of variables involved in the system as 
follows. We swap d and u. The schematic of the system is now shown in 
Figure Notice that the mathematical model of the system is still given by 
Now take the controller C given by fi48\ ). Clearly, the features of the 
interconnection are still there. What the controller now does is shut 
down the pump. This controller does not keep the water level constant. But, 
that is not the fact that we are interested in. The interesting observation is 
that now the interconnection does make sense. 

These two examples suggest the following facts. 

• We may need to introduce a constraint for systems interconnection to 
make sense. The constraint cannot be formulated based on the math- 
ematical representation of the systems alone. The interconnections 
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described in Example [19] and Example [20] share the same mathemati- 
cal representation, yet in one case the constraint is not satisfied, while 
in the other it is. This is in contrast with the regularity constraint, 
where the constraint can actually be derived from the behaviors them- 
selves. 

• The new constraint is different from the regularity constraint. Exam- 
ple [19] describes an interconnection where the regularity constraint is 
satisfied, while the new constraint that we are going to formulate is 
not satisfied. 

As indicated by the Example [TUl the constraint is violated when the 
plant is restricted through a variable that is inherently an output of the 
system. That is, the variable is physically dictated to be an output of the 
system. The information that a variable is an output cannot be deduced 
from the mathematical description of the system, rather it has to be provided 
in addition to the description of the plant. System variables that have to be 
output variables by physical consideration, are called declared outputs. We 
then require that the controller accepts the declared output of the plant as its 
input, for the interconnection to make sense. To say it differently, suppose 
that y is a (set of) variable(s) that is a part of the control variables. If y 
is declared as output because of some physical interpretation of the system, 
we want to input-output partition the variables of the controller, such that 
y belongs to the input part. Input-output partitioning of the variables of a 
linear system has been introduced in Section [2] 

The control problem with input-output partitioning constraint for linear 
systems is then formally defined as follows. 

Control with input-output partition constraint. Given a control prob- 
lem, where the plant is 



The control variables are u and y, where y is the declared output variables 
of the plant. The to-be-controlled variable is w. The desired specification 
is given as 



V 



(w,u,y) | R 



Find a regular controller C described as 




(51) 



C = {(«,y)|C 1 (|)« + ft(|) y = 



(52) 
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such that C implements S and the variables in C can be input-output parti- 
tioned such that y belongs to the input part. 

5.2 The solution 

We shall now devise an algorithm that solves the problem. We assume that 
the specification S is regularly achievable (otherwise the problem is clearly 
not solvable) 

Notation 21 We denote the class of regular controllers that implements S 
as 4 cg . 

To find a solution to the problem, we need to use the following result. 
Lemma 22 Given a controller 

c = {("-^i c '(s)» +C2 (l) ! '= }- (53 » 

Without loss of generality we assume that [C\ C2] is full row rank. The 
following statements are equivalent. 

(i) The variables in C can be partitioned such that y belongs to the input 
part 

(ii) C\ is full row rank. 

(Hi) For any y G £°°(M,R y ) there exists aue £°°(R,R U ) such that (u,y) G 
C. 

Proof, (ii => i) Suppose that C\ is full row rank. If C\ is a square matrix, 
then we already have an input-output partition with u as the output and y 
as the input. If C\ is not square, then we can partition it into 

Ci = [ C n Ci a ] , (54) 

possibly after rearranging the columns, such that Cn is a square matrix 
with full row rank. We can also partition u accordingly into ui and U2. 
Now we have an input-output partition with ui as the output and U2 and 
y as the input. 

(i iii) Suppose that the variables in C can be partitioned such that y 
belongs to the input partition. This means we can partition u into ui and 
U2, such that we have ui as the output and U2 and y as the input. So we 
can partition C accordingly such that (|54|) holds. Following the elimination 
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procedure in Section [21 we can eliminate ui and find that the behavior in 
terms of y and u 2 is £°°(M,]Ry +U2 ). 

(iii ii) We shall prove it by contradiction. Suppose that C\ is not 
full row rank. The matrix [C\ C2] can be transformed (by premultiplication 
with a suitable unimodular matrix) into 

°1 °21 

.0 c 22 y 

where C[ and C 22 are full row rank. Following the elimination procedure 
in Section [21 we can eliminate u and find that the behavior in terms of y 
is the kernel of C^Jr). Hence, we cannot choose any y G <£°°(lR,]R y ) as a 
trajectory of y. ■ 

It is straightforward to see that problem can be reformulated as follows. 
Problem. Find a controller C G £^ g in the form of 

c = { ( * lS,|Cl (s)" +C2 (s) y=0 }' (55) 

where C\ is full row rank. 

We shall use the following result. 

Lemma 23 Let X be a subset of such that for any C G£^ g there exists 
a G X such that C C . Then there exists a C G £^ g that solves the 
control problem with input-output partitioning constraint if and only if there 
exists a G X that does so. 

Proof, (if) Trivial, since X C ^ eg . 

(only if) Suppose that C G £^ g satisfies the constraint. We shall show 
that any C G (£^ g such that C C also satisfies the constraint. Let C be 
given as the kernel of [ C\ C2 1 as in (|55p . We know that C\ is full row 
rank. Since C C 0, there must be a full row rank matrix F such that 
is the kernel of [ FC\ FC2 ] ■ We also know that FC\ is full row rank. 
Therefore also solves the problem. ■ 

This lemma tells us that if we can construct a subset of (£^ g with the 
property of X, we do not need to search for the candidate controller in the 
whole <~5 g - Rather, we can restrict our attention in X. Theorem [9] shows 
that £ca! has the desired property. Thus, we shall try to construct the 
desired controller in £can> which we can parametrize according to Theorem 

El 

A solution to the control problem can be found by executing the following 
steps. 
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Step 1. Construct the canonical controller C can for the problem. Since S 
is regularly implementable, we know that the canonical controller im- 
plements S. 

Step 2. Construct a controller C G Ccan- The proof of Theorem [8] contains 
information on how to construct C from a regular controller. Denote 
the kernel representation of C and the control manifest behavior, V c , 

by 



C 



d 



dt 



(«,y) I 



(56) 
(57) 



respectively. 



Step 3. Following Theorem [TT1 any controller C in (teal can be represented 
as 



C' = \(u,y) | {d + VPC 



u + (C 2 + VP 2 ) 



y = o 



The kernel representation of a controller in £ C an that satisfies the input- 
output partitioning constraint can be found by finding a matrix V such 
that C\ + VP\ is full row rank. 

The necessary and sufficient condition for the existence of such a matrix 
V is given in the following lemma. 

Lemma 24 Given polynomial matrices C G R cxq [£] and P G M pXq [£]. 
There exists a polynomial matrix V G M CXp [£] such that C + VP is full 
row rank if and only if 

' P 

C 



rank 



> c. 



(58) 



Proof, (only if) Consider the following relation 



I 




' p ' 




p 


V I 




c 




C + VP 



(59) 



Suppose that C + VP is full row rank. This means it has a rank of c. Since 
the left multiplication in (|59p is a unimodular transformation, we know that 



rank 



P 

C 



rank 



P 

C + VP 



> c. 



(60) 
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(if) Assume that (|58p holds. If C or P is zero, we can obviously choose a 
V such that C+VP full row rank. We exclude these trivial cases and suppose 
that both P and C are nonzero. Since the rank of a polynomial matrix is 
not affected by left and right multiplication by unimodular matrices, we can 
assume without any loss of generality that P has the the Smith form, that 
is, the form of 

" Pi o 




P 



(61) 



where P\ is a diagonal matrix with nonzero determinant. Furthermore, 
with some appropriate left multiplication with unimodular matrix, we can 
transform C to the following form. 



C 



Cii 


C\2 


C21 












(62) 



where C12 and C21 are full row rank. Denote the rank of P\, C12, and C21 
as p', c', and c" respectively. We have the following relation 



rank 



P 

C 



rank C\i + rank 

= c' + p', 
> c. 

Thus 

p' > c - c'. 
We can partition V accordingly to form 



Pi 

C21 



V 



V u V 12 
V21 V22 
V31 V32 



(63) 

(64) 
(65) 

(66) 



(67) 



We structure V to have the following form. 



V 








(68) 



where V31 is to be chosen later. Therefore 



C + VP 



C11 


C\2 


C21 





V 3 iPi 






(69) 
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Our goal is to make C + VP a full row rank matrix. Therefore, 



C21 
V31P1 

has to be a full row rank matrix. Since C21 is full row rank and has the 
rank of c", we can find c" columns of C21 that form a square matrix with 
nonzero determinant. Denote this selection as N, where N C {1, 2, • • • , p'}. 
We construct V31 G IR( C ~ C ~ c ) Xm [£] such that the entries on the i— th column 
of V31 are zero if i G N. The remaining (m' — c") columns of V31 form a 
(c — c' — c") by (m' — c") matrix. From (|66p we know that it is a wide 
matrix. We choose the values of the entries of these columns such that this 



wide matrix is full row rank. It follows that 



C21 
^31^1 



is a full row rank 



matrix and hence C + VP is full row rank. ■ 

To conclude, the following is the algorithm to solve the control problem 
with input-output partitioning constraint. 

Algorithm 25 The following steps provide a solution to the problem if and 
only if it is solvable. 

1. Verify if the specification S is regularly achievable. If so, go to step 2, 
otherwise the problen is not solvable. 

2. Construct the canonical controller for this problem, denote it as C can . 

3. Construct a regular controller C € Ccan • Theorem [3 guarantees that this 
can be done. The controller C and the control manifest behavior V c can be 
represented in the form of 

c = I oi (s) - + (s) » " °} ■ < 70 > 
"c = {(«,,)in(!)» + p 2 (!> = oj. in) 



4. Verify if 



rank 



Mi 
Pi 



>P(C), 



(72) 



where p(C) denotes the number of output variables of C. If (TTSUis satisfied, 
go to step 5, otherwise the problem is not solvable. 

5. Compute a V such that C\ + VP\ is full row rank. The existence of such 
V is guaranteed by Lemma\21 h , 



A controller that solves the control problem 



is given by 



C = {(«,») I [c 1 + vp 1 c 2 + yp 2 ](£)^]=o}. 



(73) 
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6 Concluding remarks 



We discuss a result in the field of behavioral control theory for linear systems. 
The main result of the paper is a parametrization of all regular controllers 
that are equivalent to the canonical controller <£ C an- This class of controllers 
has two nice properties: 

(i) All its members are regular controllers, and 

(ii) it acts as an upperbound to other regular controllers. This means, any 
regular controller is contained in an element of Cell. 

The special properties of the class £can and its parametrization is used 
to solve two control problems in the behavioral framework. The first control 
problem is related to designing a regular controller that uses as few con- 
trol variable as possible. The second problem is about designing a regular 
controller that satisfies a predefined input-output partitioning. 

The use of the parametrization of Ccai is not necessarily limited to the 
above mentioned problems. An interesting problem is, for example, to use 
the parametrization to construct a regular controller with as small MacMil- 
lan degree as possible [Tj. Such a result can potentially lead to the solution 
to the long standing problem of regular feedback implement ability |16j . 
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